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Abstract 

Mitogen-activated protein kinase (MAPK) cascades have crucial roles in the regulation of plant development and in plant 
responses to stress. Plant recognition of pathogen-associated molecular patterns or pathogen-derived effector proteins has 
been shown to trigger activation of several MAPKs. This then controls defence responses, including synthesis and/or 
signalling of defence hormones and activation of defence related genes. The MAPK cascade genes are highly complex and 
interconnected, and thus the precise signalling mechanisms in specific plant-pathogen interactions are still not known. 
Here we investigated the MAPK signalling network involved in immune responses of potato (Solanum tuberosum L.) to 
Potato virus Y, an important potato pathogen worldwide. Sequence analysis was performed to identify the complete MAPK 
kinase (MKK) family in potato, and to identify those regulated in the hypersensitive resistance response to Potato virus Y 
infection. Arabidopsis has 10 MKK family members, of which we identified five in potato and tomato (Solanum lycopersicum 
L), and eight in Nicotiana benthamiana. Among these, SX.MKK6 is the most strongly regulated gene in response to Potato 
virus Y. The salicylic acid treatment revealed that StMKK6 is regulated by the hormone that is in agreement with the salicylic 
acid-regulated domains found in the StMKK6 promoter. The involvement of StMKK6 in potato defence response was 
confirmed by localisation studies, where StMKK6 accumulated strongly only in Pofafo-v/rus-V-infected plants, and 
predominantly in the cell nucleus. Using a yeast two-hybrid method, we identified three StMKK6 targets downstream in the 
MAPK cascade: StMAPK4_2, StMAPK6 and StMAPKI 3. These data together provide further insight into the StMKK6 signalling 
module and its involvement in plant defence. 
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Introduction 

Mitogen-activated protein kinase (MAPK) cascades are con- 
served signalling modules in eukaryotes that transduce extracel- 
lular stimuli downstream from the receptors, thus mediating the 
intracellular responses. The plant MAPK cascades have pivotal 
roles in the regulation of plant development and in responses to a 
variety of stress stimuli, including pathogen infection, wounding, 
temperature, drought, salinity, osmolarity, UV irradiation, ozone 
and reactive oxygen species [1]. 

In a general model of the MAPK signalling cascade, activation 
of plasma membrane receptors activates the MAPK kinase kinases 
(MKKKs). These are serine or threonine kinases that phosphor- 
ylate a conserved S/T-X3 — 5-S/T motif of the downstream 
MAPK kinases (MKKs), which, in turn, phosphorylate MAPKs on 
threonine and tyrosine residues in a conserved T-X-Y motif of 
their activation loop [2] . Following this MKK alteration of the 
phosphorylation-dependent properties of their target proteins, the 
activated MAPKs translate the information further, which 



eventually leads to changes in, e.g., gene expression, cellular 
redox state, or cell integrity [3]. 

The genome of the model plant Arabidopsis thaliana encodes 60 
MKKKs, 10 MKKs and 20 MAPKs [4]. This indicates that the 
MAPK cascade might not simply consist of a single MKKK, 
MKK and MAPK connected together, but has the potential to be 
organised into many thousands of distinct MKKK-MKK-MAPK 
combinations, with some level of redundancy. To minimise 
unwanted cross-talk, the spatial and temporal activities of the 
different components must be strictly regulated [5]. 

Despite the potential multiplicity of MAPK cascades, only a 
small number of MAPK modules have been experimentally 
defined [6]. As the MKK family consists of a relatively small 
number of genes, their activity in different MAPK modules is 
widely dispersed [5]. In Arabidopsis, MKKs can be divided into 
four different groups (A— D) based on their sequence similarities 
[4,7]. Group A includes Arabidopsis thaliana AtMKKl, AtMKK2 
and AtMKK6. AtMKKl and AtMKK2 act upstream of 
AtMAPK4 in response to cold, salinity and pathogens [8,9]. 
AtMKK6 is involved in cytokinesis control and cell-cycle 
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regulation [10]. The group B MKKs includes AtMKK3, which 
functions upstream of AtMAPK6 in the regulation of jasmonic 
acid (J A) signal transduction [11] and is involved in pathogen 
defence [12]; overexpression of AtMKK3 leads to enhanced 
tolerance to salt and increased sensitivity to abscisic acid [13]. The 
group C MKKs include AtMKKf and AtMKK5, which act 
upstream of AtMAPK3 and AtMAPKG in the regulation of plant 
development and defence responses [1,14—16]. The group D 
MKKs include AtMKK7, AtMKK8, AtMKK9 and AtMKKlO. 
AtMKK9 is involved in ethylene signalling [17] and in leaf 
senescence [18], while AtMKKJ is involved in plant basal and 
systemic resistance [19]. 

During pathogen attack, MAPK signalling is an indispensable 
component of the host defence response, in a way that it is 
involved in the crosstalk between secondary messengers and 
hormones [20]. The key hormones in plant biotic interactions 
include the salicylates, jasmonates and ethylene [21], whereby 
their specific roles depend on the particular host-pathogen 
interaction. Many studies have indicated that salicylic acid (SA) 
is a key regulatory compound in disease resistance against fungi, 
bacteria and viruses (reviewed in [22]). The importance of SA in 
viral multiplication and symptom development has also been 
confirmed in potato (Solanum tuberosum L.) - Potato virus Y (PVY) 
interaction [23,24]. Depending on the virus, SA can induce 
inhibition of viral replication and cell-to-cell or long-distance viral 
movement (reviewed in [25]) and in agreement with this, SA is also 
a key component in the directing of events during and following 
hypersensitive resistance (HR) [21,24]. Hypersensitive resistance is 
an efficient defence strategy in plants, as it restricts pathogen 
growth and can be activated during host, as well as non-host, 
interactions. It involves programmed cell death and manifests in 
necrotic lesions at the site of pathogen attack (reviewed in [26]). 
Potato virus Y is a member of the Potyviridae family and it is an 
important potato pathogen worldwide. In potato, HR is conferred 
by the Ny genes (reviewed in [27]). The potato cultivar (cv.) Rywal 
carries the Ny-1 gene and it develops HR that is manifested as 
necrotic lesions on leaves 3 days following their inoculation with 
various PVY strains [28] . 

To date, a large number of members of the MAPK cascades 
from different species have been investigated, although to the best 
of our knowledge, there has been no systematic investigation of the 
MKK family and its function in defence signalling in potato. 
Moreover, no MAPK immune response network module has been 
defined for potato - PVY interactions. 

We thus first performed sequence analysis of the complete 
MKK gene family in potato, and of its close relatives, where their 
genomes have been sequenced. Based on the present transcrip- 
tome data, StMKK6 was identified as the most responsive 
member after viral attack. We further investigated the role of 
StMKK6 in the response to PVY infection at the gene expression 
level, and studied its intracellular localisation and identified its 
downstream targets in the MAPK signalling cascade. 

Materials and Methods 

Bioinformatics analysis 

Arabidopsis thaliana MKK gene family (The Arabidopsis 
Information Resource; TAIR: http://www.arabidopsis.org/) was 
blasted (tBLASTx algorithm) [29] against the potato [30], tomato 
(Solanum lycopersicum L.) and Nicotiana benthamiana (The SOL 
Genomics Network: http://solgenomics.net/) genomes to identify 
the MKK gene family in all three Solanaceae species. The names 
of the potato MKK genes were assigned based on the names of 
their apparent A thaliana orthologues. 



To identify and gather all of the available sequence information on 
StMKK6 (GenBank accession number KF837127), StMAPK4_l 
(GenBank accession number KJ027594), StMAPK4_2 (GenBank 
accession number KJ027595), StMAPK6 (GenBank accession 
number KJ027596) and StMAPK13 (GenBank accession number 
KJ027597), several database searches were performed. Sequences of 
our isolated genes were blasted (tBLASTx; [29]) against NCBI 
(National Centre for Biotechnology Information: http://www.ncbi. 
nlm.nih.gov/); UniProt (Universal Protein Resource: http://www. 
uniprot.org/); TAIR; PlantGDB (http://www.plantgdb.org/); POCI 
(Potato Oligo Chip Initiative: http:/ /pgrc-35. ipk-gatersleben.de/pls/ 
htmldb_pgrc/fi ) p = 194:1); The Gene Index Project (http:// 
compbio.dfci.harvard.edu/tgi/); and TIGR Plant Transcripts As- 
semblies (http:/ / plantta.jcvi.org/ search.shtml). 

To identify SxMKK6 orthologues we searched the UniProt 
database using the BLASTp algorithm [29]. Multiple-sequence 
alignments were performed using the MAFFT programme 
(version 7) [31]. Phylogenetic trees were constructed based on 
protein sequences using the MEGA 5 software [32], with the 
neighbour-joining method [33]. Bootstrap values were derived 
from 1 ,000 replicates, to quantify the relative support for branches 
of the inferred phylogenetic tree. To distinguish between MKK6 
orthologues and other MKKs, other AxMKK s were included in the 
first MAFFT analysis. 

Information on the targets of AtMKK6 was collected from the 
Arabidopsis Interactome Network Map [34]. The corresponding 
targets in potato were assigned by A. thaliana-potato orthologue 
and potato paralogue connections [35] . 

The location of the StMKK6 gene (PGSC0003DMG403005720) 
within the potato genome was identified using the Ensembl Plants 
portal (http://plants.ensembl.org/index.html). The promoter se- 
quences were obtained from the Soktnum_tuberosum-PGSC_ 
DM_v34_superscaffolds database. Promoter sequences were ana- 
lysed using the PlantCARE web service [36] . 

Localisation predictions of the A. thaliana and potato MKK6 
proteins were performed using PredictProtein [37] and the data on 
the whole A. thaliana MKK family using the SUBA3 web services 
[38]. 

Metadata analysis 

The expression data of the StMKK genes in developmental 
tissues and under stress conditions were collected through the use 
of the Bio-Analytic Resource for Plant Biology (http://bar. 
utoronto.ca/welcome.htm) and the Potato eFP browser [39]. 
The developmental tissues were vegetative and reproductive 
organs from greenhouse-grown plants. For each MKK, the RNA 
expression signals (FPKM values) based on the RNA sequencing of 
double monoploid S. tuberosum Group Phureja clone DM 1-3 
(DM) [40] and the genome sequence and analysis of the tuber crop 
potato [30] were collected. hxMKK6 transcription data from 
various datasets were analysed using Genevestigator [41] and only 
two-fold changes were regarded as significant. 

We collected data from previously reported microarray 
expression analyses of the potato HR cv. Rywal and the SA 
deficient NahG-Rywal infected with PVY N " Wl ([24], GEO: 
GSE46180), for all of the potato MKKs and MAPKs. Microarray 
identifiers were connected to the PGSC gene identifiers [35] and 
only the probes with less than four single nucleotide polymor- 
phisms compared to the PGSC transcripts were used. 

Co-expression analysis of microarray expression datasets was 
carried out with Biolayout Express3D 3.0 software [42], with 
ratios of PVY/mock normalised signals. The network was 
constructed using the Pearson correlation coefficient threshold of 
0.98, and the graph was clustered using the Markov clustering 
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algorithm with inflation 2.3; the other parameters were set at the 
default values. 

Gene expression analysis 

Total RNA from control and SA- treated samples was extracted 
using MagMAX-96 Total RNA Isolation Kit (Life Technologies) 
according to manufacturer's instructions. RNA concentration was 
quantified by UV absorption at 260 nm using a NanoDrop 
ND1000 spectrophotometer (Nanodrop technologies). Integrity 
and purity of the RNA samples were determined by agarose gel 
electrophoresis and OD 260/280 nm absorption ratios. To 
confirm microarray results the same RNA samples previously 
analysed in microarray experiments [24] were used. Reverse 
transcription was performed on 1 |ig of total RNA using High 
Capacity cDNA Reverse Transcription Kit (Applied Biosystems). 
The expression of StMKK6 was assayed by real-time PCR 
(qPCR). The primers and probe (MKK6_F, MKK6_R and 
MKK6_S) targeting StMKK6 were designed by Primer Express 
2.0 (Applied Biosystems) (Table SI in File SI). The qPCR was 
performed using TaqMan chemistry, with Cq values determined 
as described previously [24]. All reactions were run on a 7500 Fast 
Real-Time PCR System (Applied Biosystems) in 5 |U.l volume with 
300 nM concentration of primers and 150 nM of probe and 
performed in duplicate. Linearity (R2) and efficiency (E = 1 0 [ — 1 / 
slope]) [43] of each reaction were compared to the accepted 
values. The standard curve method was used for relative gene 
expression quantification. Target gene accumulation was normal- 
ised to two endogenous control genes, cytochrome oxidase (COX) 
[44] and elongation factor 1 (EF-1) [45]. To determine differen- 
tially expressed genes, the Student t-test was performed. 

Cloning 

The full-length sequences of StMKK6, StMAPK4_l, 
StMAPK4_2, StMAPK6 and StMAPK13 were amplified from 
potato cv. Rywal cDNA with the primers listed in Table S 1 in File 
SI. The fragments were inserted into the pJET 1.2 blunt cloning 
vector (Thermo Scientific) and sequenced (CATC Biotech). For 
the localisation studies, StMKK6 was cloned into the pENTR/D 
TOPO vector (Invitrogen), using the LR reaction, following the 
manufacturer protocol, and recombined into the binary destina- 
tion vectors pH7YWG2 (YFP) and PH7CWG2 (CFP) [46], to 
produce proteins with C-terminal YFP/CFP fusion. 

The StMKK6 promoter was amplified using the genomic 
library from potato cv. Sante as the template, which was 
constructed using GenomeWalker Universal kit (BD Biosciences 
Clontech), as described previously [47]. In the genome walking 
PCR amplifications, the Advantage 2 Polymerase Mix (Clontech) 
was used with the PCR conditions suggested by the manufacturer. 
The adaptor primer API and a nested primer AP2, that were 
provided by the manufacturer were paired with the nested reverse 
gene-specific primers MKK6_AP1 and MKK6_AP2 for amplifi- 
cation of the region upstream from the StMKK6 gene. The 
fragment was inserted into the pJET 1.2 blunt cloning vector 
(Thermo Scientific) and sequenced. The 200-bp region of the 
StMKK6 promoter (identical in all three of the investigated 
genotypes) was afterwards cloned in front of the StMKK6 gene in 
the pH7YWG2 vector, using QuikChange II XL Site-Directed 
Mutagenesis kit (Agilent Technologies), with the MAPKprom- 
YFP_F and MKK6prom-YFP_R primers. The StMKK6 promot- 
er region from potato cv. Rywal was amplified from the genomic 
DNA using the primers listed in Table S 1 in File S 1 . 

For the yeast two-hybrid analysis, the coding regions of the 
selected genes were cloned into pGBKT7 (StMKK6) or pGADT7 
(StMAPK4_l, - 4_2, - 6, -13) using QuikChange II Site-Directed 



Mutagenesis kit (Agilent Technologies). The primers used are 
listed in Table SI in File SI. 

Hormonal treatment 

Two different potato genotypes were used in this study 
including 1 non-transgenic cv. Rywal and 1 transgenic line NahG, 
transgenic plants of cv. Rywal (NahG-Rywal) [24] . Potato plants 
were grown in stem node tissue culture. Two weeks after node 
segmentation they were transferred to soil and grown as previously 
described by Baebler et al. 2009 [45]. 

SA treatment was performed spraying plants with SA-analog 
300 uM INA (98% 2,6-Dichloroisonicotinic acid, Aldrich) in 
distilled water. Control plants were sprayed with distilled water 
alone. Leaves were harvested after 24 h of treatment and 
immediately frozen in liquid nitrogen. Two biological replicates 
per treatment were analysed. 

Localisation studies 

N. benthamiana plants were grown in a growth chamber at 22/ 
20°C under a 16/8 h light/dark cycle. PVY NTN (isolate NIB- 
NTN, AJ585342) inoculation of N. benthamiana was performed 
as described by Baebler et al. 2009 [45]. Plasmids pH7YWG2 
(YFP), pH7CWG2 (CFP), P H7YWG2::StMKK6-YFP, 
P H7CWG2::StMKK6-CFP under the 35S promoter, and 
pH7YWG2::StMKK6 under the control of the native promoter, 
were introduced into Agrobacterium tumefaciens strain GV3101 
by electroporation. The A. tumefaciens cells were cultured, 
harvested by centrifugation, and re-suspended in 0.2 mM 
acetosyringone solution, and infiltrated into 4— 5-week-old N. 
benthamiana leaves, 8 days after PVY/mock inoculation. 

After 48 h to 72 h, YFP/CFP were visualized with a Leica TCS 
SP5 laser-scanning microscope mounted on a Leica DMI 6000 CS 
inverted microscope (Leica Microsystems, Germany), with an HC 
PL FLUOTAR 10.0x0.30 DRY objective. For excitation, the 
458 nm (enhanced cyan fluorescent protein) and 514 nm 
(enhanced yellow fluorescent protein) lines of an Argon laser were 
used. The enhanced cyan fluorescent protein emission was 
measured from 475 nm to 495 nm, and the enhanced yellow 
fluorescent protein from 525 nm to 550 nm. Autofluorescence 
emission was measured from 690 nm to 750 nm. Differential 
interference contrast (DIC) images were captured using the 
transmission light detector of the confocal microscope. All of the 
images were acquired with the 10 x objective or with an additional 
3.64 x zoom. The acquired images were processed using the Leica 
LAS AF Lite software (Leica Microsystems). pH7YWG2 and 
pH7CWG2 without the inserted StMKK6 were used as the 
controls. 

Yeast two-hybrid assay 

The coding regions of selected genes were cloned into pGBKT7 
{StMKKQ or pGADT7 {StMAPK4_l, -4_2, -6, -13). To 
introduce insertions into plasmids the QuikChange II Site- 
Directed Mutagenesis Kit (AgilentTechnologies) was used. The 
primers used are listed in Table SI in File SI. 

Yeast two-hybrid analysis was performed using the Matchmaker 
GAL4-based two-hybrid assay (Clontech). The vector pGBKT7, 
containing bait protein fused to Gal DNA-binding domain, was 
used to transform the yeast strain Y2H Gold through the 
polyethylene glycol/LiAc-based method (Clontech), according to 
manufacturer instructions. After selection on SDALeu medium, 
the transformed colonies were used for co-transformation follow- 
ing the protocol described by Clontech, with the pGADT7 vector. 
It was used as a prey and contained target genes fused to Gal 
DNA-activation domain. Co-transformants were selected on SD/- 
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Leu/ -Tip (DDO) medium. Positive interaction transformants 
were selected on SD/-Leu/-Trp/-His/-Ade/x-a-Gal/Aba 
(QDO/X/A) medium. All of the genes were previously tested 
for autoactivation and toxicity, plating yeast colonies transformed 
with pGBKT7_StMKK6 on SD/-Trp/x-a-Gal/Aba and 
pGADT7_StMAPK4_l, -4_2, -6, -13 on SD/-Leu/x-a-Gal/ 
Aba. Positive (pGBKT7_SV40 large T-antigen -pGADT7_53) 
and negative (pGBKT7_SV40 large T-antigen -pGADT7_Lam) 
interaction controls, included in the Matchmaker Gold Yeast 
Two-Hybrid System (Clontech), were performed in parallel. 

Results 

The MKK family in potato 

Our genome analysis of the potato, tomato and Nicotiana 
attenuata genomes identified five MKK genes in each, with four in 
the tobacco Nicotiana tabacum L. genome, and eight in the N. 
benthamiana genome (Figure 1). Based on the phylogenetic tree of 
Arabidopsis and five Solanaceae species, for only two of the A. 
thaliana genes (MMKK3, AtMKK6) a specific orthologue could be 
assigned from potato, tomato or Nicotiana spp., while the other 
MKKs form larger orthologue groups. The MKK family in potato, 
tomato and Nicotiana spp. can also be divided into four groups: 
group A (A. thaliana and N. benthamiana, three genes; potato, 
tomato, N. tabacum and N. attenuata, two genes); B (N. benthamiana, 
two genes; the rest, one gene each); C (A. thaliana and N. 
benthamiana, two genes; potato, tomato, N. tabacum and N. 
attenuata, one gene each); and D (A. thaliana, four genes; Solanaceae 
spp., one gene each, except JV. tabacum which has none). Out of the 
five examined Solanaceae spp., N. benthamiana is the only one of 
these in which duplication of the MKK genes occurred after 
speciation and thus N. benthamiana has possible paralogues for 
MKK3 (NbS00014007gO022.1, NbS00037566g0012.1), MKK4/5 
(NbS00012713g0030.1, NbS00006609g0002.1) and MKK6 
(NbS00006857g0015.1, NbS00045036g0007. 1) (Figure 1). 

Regulation of potato MKK family gene expression in 
developmental processes and stress 

We first analysed the gene-expression metadata for the potato 
MKK family for developmental processes and after exposure to 
different biotic and abiotic stress [40]. The tissue specificities of the 
MKK family members differed between both of the potato 
varieties studied here, the potato Group Phureja Clone DM, and 
Group Tuberosum Clone RH, especially regarding the expression 
in tubers and stolons. In general, this analysis showed that in 
potato, StMKKl/2 and StMKK4/5 are relatively uniformly 
expressed, with the highest expression in tubers and stolons. 
StMKK3 transcripts had the highest abundance in carpels and 
petals, and StMKK7/9 in tubers and flowers, and also in roots for 
Clone RH. Under biotic and abiotic stress conditions, all of the 
potato MKKs except StMKK4/5, were down-regulated after 
wounding. None of the other stress conditions and hormonal 
treatments influenced the MKK family gene expression by more 
than two-fold compared to controls (File S2). 

In a previous transcriptome analysis of potato HR responses to 
PVY infection [24], we checked the expression profile of identified 
potato MKKs [35]. According to the gene expression profiles 
obtained for the potato cv. Rywal, three out of five MKKs showed 
differential expression for the HR response to PVY: StMKK3, 
StMKK6 and SMKK4/5 (Figure 2). StMKK6 showed the 
strongest response of all, thus its expression data were also 
confirmed by quantitative real-time PGR (qPCR) (Table S2 in File 
SI). StMKK6 was strongly up-regulated before lesion formation (3 
dpi). In contrast, at the same time, StMKK3 was down-regulated. 



StMKK4/5 showed two-fold induction at the later time (6 dpi) 
following inoculation. For the SA-deficient plants (NahG-Rywal), 
the regulation of these MKKs at the gene level was either 
attenuated (for StMKfO and StMKK4/5) or delayed (for 
StMKK6). Only StMKK7/9 was induced earlier and stronger in 
the NahG-Rywal plants, compared to the non-transgenic plants. 

Based on these data we focused on further investigations on the 
AtMKK6 orthologue in potato, StMKKS. 

StMKK6 orthologues across the plant kingdom 

To study the sequence diversity and evolution of the MKK6 
gene, we performed a phylogenetic analysis within the plant 
kingdom. We identified 39 plant MKK6 orthologues among 
higher plants. As well as for different families of Angiosperms, one 
MKKS gene was also identified in Gymnosperms, and one in 
mosses (Figure 3 and File S3). MKK6 is mostly a single-copy gene, 
as these 39 MKK6 orthologues belong to 31 plant species 
(Figure 3). 

Comparison of the MKK6 sequence from potato cv. Rywal and 
the published potato genome sequence of S. tuberosum Group 
Phureja, Clone DM showed that the coding domain sequences are 
different for three nucleotides at different sites, which results in 
only one amino-acid difference; this is not part of the protein 
kinase domain. 

Regulation of StMKK6 gene expression in development 
and stress 

To better understand the function of the MKK6s, we examined 
the gene expression data that is available for A. thaliana MKK6 in 
Genevestigator. AtMKK6 is up-regulated in developmental 
processes (callus formation, germination) and after treatments 
with translation inhibitor cyclohexamide and hormones auxins. 
AtMKK6 is also up-regulated after flagellin 22 treatment, and 
under abiotic stress (cold, salt, acidic pH, hypoxia). AtMKK6 is 
down-regulated after treatments with abscisic acid and combina- 
tions of abscisic acid with jasmonate or SA. 

The expression profiles of StMKK6 obtained from the Potato 
eFP Browser database showed that in the potato Group Phureja 
Clone DM, StMKK6 tissue expression was high in the tubers and 
callus, and very low (below detection levels) in the flowers. In the 
Group Tuberosum Clone RH under non-stress conditions, 
StMKK6 gene expression was strong also in the shoot apex 
(Figure 4A; [40]). Among all of the stress-related conditions 
included, compared to the controls, StMKK6 was down-regulated 
after wounding of the secondary leaves and after treatment with 
the fungal elicitor DL-(3-amino-n-butyric acid (Figure 4B). 

Analysis of the SX.MKK6 native promoter 

In the Potato Genome Sequencing Consortium (PGSC) Potato 
Genome Browser, StMKK6 is assigned to PGSC gene ID 
PGSC0003DMG403005720, and it is localised to chromosome 
3. Based on this information, we also cloned the StMKK6 
promoter sequence from two cultivars of S. tuberosum Group 
Tuberosum, both of which are resistant to PVY infection: cv. 
Rywal (GenBank accession number KF837128) and cv. Sante 
(GenBank accession number KF837129). The 790 bp-long 
promoter isolated from cv. Rywal (Table S3 in File SI) is identical 
to the promoter of S. tuberosum Group Phureja, but different from 
that of cv. Sante (Table S4 in File SI), due to a 27-bp-long region 
that is localised 230 bp upstream ofStMKK6. In general, the most 
common domains were the TATA-box, CAAT box and domains 
related to light responsiveness. Based on the development and 
stress-related domains found, StMKK6 appears to be involved in 
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Figure 1 . Phylogenetic tree of /WATf family in A. thaliana and five Solanaceae species. The species are A. thaliana (At), potato (Sotub), tomato 
(Solyc), N. benthamiana (Nb), N. attenuata (Na) and N. tabacum (Nt). Genes are grouped into 4 groups: A (green), B (red), C (blue) and D (orange) [4]. 
Potato genes are marked with dots. The numbers on the nodes are percentages from a bootstrap analysis of 1000 replicates. The scale bar indicates 
the branch length that corresponds to 0.06 substitutions per site. 
doi:10.1371/journal.pone.0104553.g001 



development of the endosperm and in responses to heat stress and 
wounding. StMKK6 is also predicted to be involved in gibberellic 
acid and salicylic acid signalling, and is under regulation by the 
circadian clock (Figure 5). 

Regulation of StMKK6 by salicylic acid 

Promoter analysis of StMKK6 has demonstrated potential 
involvement of StMKK6 in SA signalling. Therefore, we 
examined expression of StMKK6 under the SA treatment. The 
expression of StMKK6 was measured by qPCR 24 h after 
treatment with SA in Rywal and NahG-Rywal plants. 

The treatment showed StMKK6 to be sensitive to the present 
active form of SA since significant changes of StMKK6 expression 
were observed in the SA-deficient genotype NahG-Rywal 
(Figure 6). There was a substantial, two-fold rise, of the StMKK6 
RNA, compared to the non-treated plants. In the Rywal plants, 



that received the same SA treatment, the change in the StMKK6 
expression was not significant (Table S5 in File SI). 

We also observed significant difference in the RNA level of the 
endogenous StMKK6 in the non-treated Rywal, compared to the 
NahG-Rywal (Figure 6). In Rywal plants there was at least twice 
as much of expressed StMKK6 comparing to the level in the 
NahG-Rywal (Table S6 in File SI). 

Localisation of the StMKK6 protein in healthy and 
PVY-infected epidermal leaf cells 

To learn more about the function of the StMKK6 protein in 
defence responses against the virus PVY, we analysed StMKKG 
localisation in healthy and infected N. benthamiana plants. The 
infected leaves of N. benthamiana developed no necrotic lesions 
due to infection with PVY or overexpression of StMKK6. In-silico 
analysis predicted the localisation of StMKK6 as in the nucleus 
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Figure 2. Gene expression pattern of MKK family in the HR response against PVY. Cultivar Rywal (HR response, conferred by Ny-1 gene) 
and NahG-Rywal (impaired accumulation of SA) were analysed for whole transcriptome response 1, 3 and 6 days after PVY N " W| infection [24], A. 
thaliana and S. tuberosum PGSC orthologues were assigned to each probe. Log 2 fold changes of PVY in infected vs. mock-inoculated plants are 
indicated for each time point (1, 3 and 6 dpi). Statistically significant differences (FDR corrected p<0.05) are in bold. Up-regulated values are in blue 
and down-regulated values are in yellow. 
doi:1 0.1 371 /journal.pone.01 04553.g002 
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Figure 3. Phylogenetic tree of potato MKK6 and its orthologues in different plant species, divided into classes. The species of origin for 
each MKK is marked by species acronym before the protein name: Cr, Capsella rubella; Aly, Arabidopsis lyrata; At, Arabidopsis thaliana; Br, Brassica 
rapa; Bn, Brassica napus; Gh, Gossypium hirsutum; Pp, Prunus persica; Mt, Medicago trunculata; Gm, Glycine max; Pt, Populus trichocarpa; Rc, Ricinus 
communis; Vv, Vitis vinifera; Oo, Origanum onites; St, Solanum tuberosum; Le, Solanum lycopersicum (Lycopersicum esculetum); Nt, Nicotiana tabacum; 
Na, Nicotiana attenuata; Nb, Nicotiana benthamiana; Ma, Musa acuminata; Hv, Hordeum vulgare, Ata, Aegilops tauschii; Tu, Triticum urartu; Bd, 
Brachypodium distachyon; Os, Oryza sativa, Og, Oryza glaberrima; Ob, Oryza brachyantha; Si, Setaria italica; Sb, Sorghum bicolor; Zm, Zea mays, Ps, 
Picea sitchensis; Sm, Selaginella moellendorffii. Potato MKK6 is marked with red dot. The numbers on the nodes are percentages from a bootstrap 
analysis of 1000 replicates. The scale bar indicates the branch length that corresponds to 0.06 substitutions per site. 
doi:1 0.1 371 /journal.pone.01 04553.g003 



and cytoplasm. The same results were obtained also for AtMKK6 
(Table S7 in File SI). An in-vivo study of the subcellular 
localisation of StMKK6 was performed under the native promoter 
and in C -terminal translational fusion with YFP in JV. benthami- 
ana leaves. In the mock-inoculated leaves, no fluorescence of the 
StMKK6-YFP fusion protein was observed (Figure 7A), while in 
PVY-infected leaves, there was strong StMKK6-YFP fluorescence 
in the nucleus and very weak StMKK6-YFP fluorescence in the 
cytoplasm (Figure 7B). 

The fluorescent protein itself did not affect the localisation of 
StMKK6, as fusions with different fluorescent proteins were 
observed in the same intracellular compartments. The overex- 
pression of StMKKG, driven by the CaMV 35S promoter, resulted 
in localisation of the protein in the cytoplasm and nucleus in mock- 



inoculated as well as in PVY-infected leaves (Figures SIB and 
SIC). 

Search for SX.MKK6 co-regulated genes and potential 
targets in the MAPK cascade 

To further improve our knowledge of StMKK6 function, we 
searched for the other components of the MAPK module that is 
involved in the HR response in potato, the genes that are co- 
expressed with StMKK6, and the possible MAPK targets. 

We found 22 co-expressed genes (File S4), for which we 
analysed the promoter regions. The most common biotic stress- 
related domains in their promoters were TC-rich repeat motifs (an 
element involved in defence and stress responsiveness; 13 
promoters with the region), TGACG and CGTCA motifs 
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Figure 4. StMKK6 gene expression in development and stress [40]. Data are obtained from the potato eFP browser [39]. A. Tissue and 
developmental gene expression pattern of StMKK6 in double monoploid S. tuberosum Group Phureja (left) and in heterozygous diploid 5. tuberosum 
Group Tuberosum (right). B. Changes in StMKK6 gene expression under biotic stress conditions (treatments of leaves with P. infestans, SA analogues 
acibenzolar-S-methyl and fungal elicitor DL-fj-amino-n-butyric acid and wounding) and abiotic stress conditions (treatments of whole plants with salt, 
heat and hormones cytokinins, 6 benzylaminopurine; gibberellins, GA3; abscisic acid; and auxin, IAA) and control treatments. 
doi:10.1371/journal.pone.0104553.g004 



(regulatory elements involved in methyl jasmonate responsiveness; 
9 promoters with the region), Box-Wl (fungal elicitor responsive 
element; 8 promoters with the region), and TCA elements (SA 
responsiveness; 8 promoters with the region). The only region that 
was common to all 22 genes was the Skn-1 motif, which is related 
to endosperm expression. 

In an A. thaliana interactome study [34], AtMKK6 was shown 
to interact with: AtMAPK4 (AT4G01370), AtMAPK6 
(AT2G43790), AtMAPKl 1 (AT1G01560), AtMAPK12 
(AT2G46070), AtMAPKl 3 (AT1G07880), and the calmodulin- 
like protein TCH3 (AT2G41100). 

We first identified the orthologues of all of the potential 
interacting MAPKs in potato. AtMAPK4 together with At- 
MAPKl 1 and AtMAPKl 2 represent one orthologue cluster with 
StMAPK4_l (PGSC003DMP400037535) and StMAPK4_2 
(PGCS003DMP400000144). We also identified specific ortholo- 
gues of AtMAPK6 in Sotub08g0 10260. 1.1, and for AtMAPKl 3 in 
PGSC003DMP400044003 (Figure S2). We further cloned all of 



the potential MAPK interactors of StMKK6 from potato cv. 
Rywal. Comparisons of the four MAPKs cloned from potato cv. 
Rywal and those published in the potato genome sequence of 5. 
tuberosum Group Phureja, Clone DM showed Rywal MAPK4_1 
to be one-amino-acid different from PGSC0003DMP400037535, 
while Rywal MAPK 4_2 and MAPK 13 are identical to their 
Phureja orthologues. Rywal MAPK6 had no corresponding 
proteins in Group Phureja, Clone DM. Yeast two-hybrid assays 
confirmed the interactions of StMKKb with three MAPKs: 
StMAPK4_2, StMAPK6 and StMAPK13 (Figure 8). 
StMAPK4_l could not be confirmed as a target of StMKK6 in 
potato cv. Rywal. 

StMAPK4_2 and AtMAPKl 3 show differential gene expression 
patterns after PVY inoculation (Figure 9), and have different 
expression patterns than StMKK6 (Figure 2). 
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Figure 5. Predicted regulatory domains of SXMKK6 native promoter from S. tuberosum cv. Rywal, 790-bp-long. The crucial domains that 
regulate stress responses and development are CCAAT-box (MYBHvl binding site; yellow), GCN4 motif (endosperm expression; pink), HSE (heat 
stress; brown), P-box (gibberellins; blue), Skn-1 motif (endosperm expression; purple), TCA element (salicylic acid; green), WUN motif (wounding; red) 
and circadian motif (control of circadian clock; grey). The arrow indicates the predicted beginning of the 5' UTR region and the red shaded ATG 
indicates the start of StMKK6 coding region. The complete list of StMKK6 promoter domains from cv. Sante and cv. Rywal with their position, strand 
and sequence are specified in Tables S3 and S4 in File SI. 
doi:10.1371/journal.pone.0104553.g005 
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Figure 6. Expression of SXMKK6\n control and SA-treated Rywal and NahG-Rywal potato plants. Expression of the endogenous StMKK6 in 
NahG-Rywal is at least two-fold lower than in Rywal. After the SA-treatment the expression of SX.MKK6 in NahG-Rywal increased for at least three 
times, comparing to the non-treated control. Difference in the endogenous SXMKK6 expression between Rywal and NahG-Rywal as well as difference 
in the StMKK6 expression between SA-treated and non-treated NahG-Rywal plants are statistically significant (p<0.05). 
doi:10.1371/journal.pone.0104553.g006 



Discussion 

To shed further light on the hypersensitive resistance signalling 
responses during PVY infection of potato, we analysed the 
involvement of potato MKKs in this process. The results revealed 
that the StMKK6 is regulated by the SA signalling pathway in the 
HR to PVY infection. The virus induces the expression of 
StMKK6 by the protein accumulating preferentially in nucleus. 

Potato has five MKKs 

To date, MKK genes have been analysed in several plants: in A. 
thaliana, 10 have been identified [4], in rice, 8 [48], in 
Brachypodium, 12 [49], in poplar, 10 [50], in apple, 9 [51] and 
in canola, 7 [52]. In the Solanaceae family, the MKK family is 
smaller than in other plant species. Our analysis showed that the 
potato and tomato genomes have five MKK genes, and that the 
N. benthamiana genome has eight. Five MKKs were identified also 
in N. attenuata, while in N. tabacum, an additional NbMKKl 
orthologue was not found among the currently existing sequences 
(Figure 1, [53]). The phylogenetic tree of the different MKK genes 
shows that only N. benthamiana has three possible pairs of 
paralogues, with none found in potato or tomato. Higher numbers 
of MKK family members in other plant species might be the result 
of gene duplications in the MKK ancestors that occurred after 
divergence into different families. Consequently, it was not 
possible to assign a single A. thaliana orthologue to three out of 



five potato MKKs. This issue of assigning orthologues was 
previously discussed by Doczi et al. 2012 [54], and was attributed 
to the high expansion and functional diversification of the 
MAPKs, and to the rapid evolution of MAPK signalling. 

It has been speculated that all 1 0 AtMKKs are fully functional, 
although there have been doubts expressed for AtMKK8 and 
AtMKKlO. Although AtMKK8 has all of the typical MKK 
motifs, there is no evidence of its expression, while AtMKKlO 
lacks a correcdy constructed target site in the activation loop [7]. 
Thus, in plants with lower numbers of MKK genes, the plasticity 
of the MAPK signalling module is not necessarily also propor- 
tionally lower. 

StMKK6 preferentially accumulates in the nucleus only 
after PVY infection 

Transcriptomic analysis indicated that the StMKK6 gene is 
involved in the HR response against PVY, as the most intensively 
regulated of the MKKs (Figure 2). Thus, we focused on different 
aspects of StMKK6 function in this process. Here, meta-analysis of 
the publicly available datasets (Figure 4) and our in-planta 
localisation study (Figure 7) show that StMKK6 expression is very 
low in non-stressed potato leaves. Interestingly, the infection with 
PVY induces the expression of StMKK6, with the StMKK6 
protein concentrated predominandy in nucleus. 

The information on the intracellular localisation of StMKK6 
and its orthologues is currendy really limited. The only report on 
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Figure 7. Epidermal cells of N. benthamiana expressing translational fusion of StMKK6 with YFP under native promoter. Leaves were 
agroinfiltrated when the virus has spread uniformly through the inoculated leaves (8 days after inoculation) and observed after 72 h in two 
independent experiments. Examples from two plants (left and right panels) are shown. Control of transformation (fluorescent marker without StMKK6 
fusion) is in the Figure S1 A. A. Localisation of StMKK6 in mock-inoculated leaves. No fluorescence was observed. B. Localisation of StMKK6 in PVY- 
inoculated leaves, where the protein accumulates predominantly in nucleus. Additional images of StMKK6 localisation under native promoter are in 
Figures SID and SI E. 
doi:1 0.1 371 /journal.pone.01 04553.g007 



experimental results of MKKs intracellular localisation showed 
that AtMKK6 accumulated in the equatorial plane of the 
phragmoplast in dividing root epidermal cells, while all of the 
other AtMKKs (AtMKKl-5) did not [10]. In N. benthamiana, 
localisation of MKK1 (an orthologue of AtMKK7, -8 and -9) was 
studied for the HR to Phylophtora infestans infection [55,56]. N. 
benthamiana MKK1 preferentially accumulated in the nucleus 
also under non-stress conditions. Liang et al. 2013 [52], studied 
the localisation of MKK2, -3 and -4 in Brassica napus epidermal 
leaf cells, and detected the accumulation of these MKKs both in 
the cytoplasm and the nucleus. On the other hand, other reports 



on the localisation of the MKK target MAPKs show preferential 
accumulation of those proteins in the cytoplasm [52,55]. One 
however has to note that all above mentioned experiments were 
studying localisation of MKKs when expressed under strong (e.g. 
CaMV 35S) promoter, while we have studied localisation under its 
native promoter (Figure 7). Also in our hands expression of protein 
under 35S promoter caused ubiquitous cellular localisation (e.g. 
cytoplasm and nucleus, see Figures SIB and SIC) presumably due 
to excessive amounts of protein produced, which is known to 
generate localisation artefacts [57]. Anyhow, our results support 
the hypothesis that activation of MKK6 targets occurs predom- 
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Figure 8. Yeast two-hybrid assays screening StMKK6 interaction partners. StMKK6 protein was fused with Gal DNA-binding domain as bait 
and StMAPK4_1, StMAPK4_2, StMAPK6 and StMAPK13 were used as prey fused with Gal DNA-activation domain. Interaction pairs P53/SV40 large T- 
antigen and Lam/SV40 large T-antigen were used as positive and negative controls respectively. StMKK6/empty vector pair was used as a control to 
discard auto-activation of bait protein. Serial dilutions of interaction pairs were plated on DDO media for co-transformation selection. Interactions of 
bait and prey proteins were examined by assessing growth on several selective media with different levels of restrictiveness i.e. SD/-Leu/-Trp/-His 
(TDO), SD/-l_eu/-Trp/x-a-Gal/Aba (DDO/X/A) and QDO/X/A. Only co-transformed colonies growing on the most restrictive media (QDO/X/A) were 
considered as positive interaction transformants. 
doi:10.1371/journal.pone.0104553.g008 



inately in the nucleus as indicated also by in silico localisation 
prediction of identified target MAPKs and that its activitiy is 
regulated on the level of MKK6 transcription. 

Several studies have focused on the gene expression profiles of 
StMKK6 orthologues in other species. A. thaliana AXMKK6 and 
maize MKK6 (ZmMEKl) are involved in the regulation of 
cytokinesis [10,58,59] and are required for lateral root formation 
[60]. ZmMEKl is also induced by polyethylene glycol, abscisic 
acid and SA, and it is down-regulated by NaCl [61]. In rice, 
OsMEKl (referred to as OsMKK6 in Hamel et al. 2006 [7]) was 
reported to be involved in signalling of moderately low temper- 
ature stress [62,63]. The present investigation into StMKK6 
expression in different tissues of the S. tuberosum group Phureja 
shows that StMKK6 is highly expressed in tissues with intensive 
cell proliferation, similar to what has been reported for its 
orthologues in A thaliana, tobacco and maize [10,58,59,64]. 

After the SA treatment of the NahG-Rywal, the level of 
StMKK6 transcripts three-fold increased, while in non-transgenic 
Rywal we detected trend of decrease in transcription of StMKK6 
gene (Figure 6). Higher basal SA levels in cv. Rywal could explain 
different regulation of StMKK 6 after exogenous treatment. In fact, 
our results showed that basal levels of StMKK6 transcripts in 
Rywal are at least two-fold higher than in tha SA-deficient NahG- 
Rywal which indicates that SA is necessary for the maintenance of 
basal levels of the StMKK6. This effect was previously reported for 
other SA-dependent defence-related genes in cv. NahG-Rywal 
[24] and cv. NahG-Desiree [23]. The network of SA signalling is 
complex [23,65-67] and expected to show some nonlinear 
behaviour. Consequently the response we see might be both time 
as well as concentration dependent and we will not be able to fully 



predict the outcome of SA signalling before its appropriate 
dynamic model exisits. In line with hypothesis that StMKK6 is 
regulated by SA, the analysis of the StMKK6 promoter revealed 
motifs that are responsive to SA (Figure 5). To the best of our 
knowledge, there has been only one study linking StMKK6 
expression and SA signalling for any of the StMKK6 orthologues 
{ZmMEKl) [61]. In B. napus, MKK1, MKK2, MKK4 and 
MKK9 were induced by SA, while only MKK3 was not [52], 
which indicates that the majority of MKKs are targets of the SA 
signalling network. 

The microarray results show that only MKK3 (group B MKKs) 
responded in the opposite direction; i.e., MKK3 was repressed 
after viral infection, contrary to all of the other potato MKKs, 
which were induced. Similarly, it has been shown also for A. 
thaliana and tobacco, that the MKK3 acts as a negative regulator 
in the JA signalling pathway and response to herbivores 
[11,16,68]. 

MAPK signalling network in hypersensitive resistance 
response to PVY 

We identified here three downstream targets of StMKTC6 in 
potato, StMAPK4_2 (orthologue of AtMAPK4, -11 and -12), 
StMAPK6 (orthologue of AtMAPK6 and tobacco SIPK) and 
StMAPK13 (orthologue of AtMAPK13 and tobacco NTF6/ 
NRK1). All of these interactions had already been confirmed in A. 
thaliana [34,69,70]. In rice, OsMEKl (or OsMKK6) interacts 
with OsMAPKl, -3, -5 and -6 (orthologues of AtMAPK6, -3 and - 
11/4, respectively [48,63]). To date, only one MAPK interactor 
has been identified for the MKK6 orthologues in Solanaceae: in 
N. tabacum, NQKl/NtMEKl was reported to interact with 
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Figure 9. Gene expression pattern of confirmed StMKK6 targets in the HR response against PVY. The experimental setup is as in 
Figure 2. To each probe an A. thaliana orthologue, potato orthologue and a PGSC gene ID were assigned. Log2 of gene expression differences 
between PVY-infected and mock-inoculated plants are indicated for each time point (1, 3 and 6 dpi). Statistically significant differences (FDR corrected 
p-value <0.05) are given in bold. The star (*) for StMAPK6 PGSC gene ID indicates that the gene was not predicted in PGSC gene model. Up-regulated 
values are in blue and down-regulated values are in yellow. 
doi:10.1371/journal.pone.0104553.g009 



PLOS ONE | www.plosone.org 



11 



August 2014 | Volume 9 | Issue 8 | e104553 



Potato (Solanum tuberosum L) MKK6 in Response to Potato virus Y 



NTF6/NRK1 (orthologue of AtMAPK13 [64]). According to this, 
we can conclude that the MKK6 signalling module is evolution- 
arily relatively stable, as most of the interactions appear to be 
conserved between unrelated species. There are, however, some 
specificities. For example, out of two potato MAPK4/ 11/12 
paralogues StMKK6 interacts with only one of them, while in A . 
thaliana it interacts with all three, and in rice, with two. 

Among the StMKK6 downstream targets, the present study 
shows that only StMAPK13 expression (Figure 9) is significandy 
regulated in the HR response to PVY infection, although while 
StMKK6 is strongly up-regulated, StMAPK13 is strongly down- 
regulated. Interestingly StMAPK4_2 and StMAPK13 are signif- 
icantly induced in SA deficient plants at the same time points as 
StMKK6. All of the target MAPKs are, however, expressed at 
(substantially) higher levels in both the mock-infected and PVY- 
infected leaves, compared to StMKK6 (Figure S3). Therefore, we 
can hypothesise that the crucial regulation of StMAPK4_2, 6 and 
13 is not at the gene expression level, but at a later step; e.g., 
phosphorylation by newly synthesised StMKK6, or any other post- 
translational process (e.g. protein translocation or degradation). 

As well as StMKK6, we showed that StMKK4/5 (also named as 
StMEKl in some studies) is also up-regulated in the HR response 
of potato to PVY, albeit to a lower extent (slightly below the strict 
significance cut-off), and showed a delay comparable to StMKK6 
in the NahG plants (Figure 2). In N. tabacum, both NtMEK2 
(orthologue of SXMKK4/5) and NQK1 /NtMEKl (orthologue of 
StMKK6) are required for N-mediated resistance against tobacco 
mosaic virus [71-73]. Similarly, the silencing of AtMEKl and 
AtMEK2 allowed for higher amplification of Cucumber mosaic 
virus mA. thaliana [74]. This indicates the close connectivity and 
interdependency between the StMEKl and SxMEK2 signalling 
modules in plant defence. 

Since an effective HR takes place in cv. Rywal in response to 
PVY, it is possible that some of the regulated potato MKKs are 
part of this response. A comprehensive analysis of the different 
MKK functions was performed in Pto-mediated resistance in 
tomato by Ekengren et al., 2003 [75], and Pedley and Martin, 
2004 [76] . Ekengren et al. [75] showed that silencing of either 
MEK1 or MEK2 breaks the resistance against Pseudomonas 
syringae. However, Pedley and Martin [76] showed that LeMKK2 
(orthologue of AxMKK4/5) and LeMKK4 (orthologue of 
AtMKKK7/8/9) caused programmed cell death when overex- 
pressed, while LeMKKl (orthologue of AtMKKl/2) and 
LeMKK3 (orthologue of MKK6) did not. In A. thaliana, two 
studies have shown that as well as the MKK4/MKK5 (MEK2) 
module, another branch of the group A MAPKKs, AtMKKl and 
AtMKK2, function together as the second MAPK module 
involved in the HR response [68,77]. In N. benthamiana MKK1 
(orthologue of AtMKK7/8/9) is a potent inducer of HR-like cell 
death in response to P. infestans [11,55,78] and the same was 
shown also for A\MKK7 and MKK9 in A. thaliana [79]. In A. 
thaliana, MKK5 was shown to have a role in the cascade that 
triggers the HR response [80]. In the present study, there was 
strong repression of potato MKK3, while this dysregulation was 
diminished in the NahG-Rywal plants (Figure 2). 

The interconnectivity between the NbMEKl (orthologue of 
AtMKK6) and NbMEK2 (orthologue of AtMKK4/5) kinase 
signalling modules was studied mechanistically by del Pozo et al., 
2004 [81]. They showed that the silencing of NbMEKl abolished 
the cell death caused by constitutively active NbMEK2. In N. 
attenuata, MEK1 (orthologue of SXMKK6) and SIPKK (ortholo- 
gue of AtMKKl/2) are involved in the regulation of the 
accumulation of 12-oxo-phytodienoic acid and JA [53], which 
from another point, supports the role of SA in MKK6 signalling, 



due to the known antagonistic cross-talk between JA and SA in 
defence pathways [82]. 

Conclusions 

Our results show that StMKK6 is an important player in potato 
HR response against PVY infection, as shown on the gene- 
expression level and protein localisation. We identified potential 
StMKK6 downstream targets and have shown that albeit the 
signalling network seems to be evolutionary relatively stable the 
fine-tuned interdependency within the broader MAPK signalling 
network might be different in different plants, and when a plant is 
exposed to different pathogens. 

Supporting Information 

Figure SI Localisation of StMKK6, under 35S promoter 
and native promoter, in epidermal cells of N. benthami- 
ana. A. Control of transformation. Epidermal cells, transformed 
with plasmids containing 35S::pH7CWG2-CFP (left) and 
35S::pH7YWG2-YFP (right) fusion. The fluorescence of the 
CFP or YFP alone (without the fusion with StMKK6) is observed 
only in cytoplasm. B. Localisation of StMKK6 fused with CFP 
(upper panel) or YFP (lower panel) with expression under the 
CaMV 35S promoter in mock-inoculated epidermal cells. The 
protein is localised in cytoplasm and nucleus. C. Localisation of 
StMKK6 fused with CFP (upper panel) or YFP (lower panel) with 
expression under the CaMV 35S promoter in PVY-inoculated 
epidermal cells. The protein is localised in cytoplasm and nucleus. 
D. Localisation of StMKK6 fused with YFP with expression under 
native promoter in mock-inoculated epidermal cells, where no 
fluorescence is observed. E. Localisation of StMKK6 fused with 
YFP with expression under native promoter in PVY-inoculated 
epidermal cells, where the protein accumulates predominantly in 
nucleus. 
(TIF) 

Figure S2 Phylogenetic tree of potato and A. thaliana 
MAPKs from group A and group B. Potato (Sotub, PGSC 
and St) and A. thaliana (At) MAPKs from group A (MAPK3, 6 
and 10) and group B (MAPK4, 5, 11, 12 and 13) as was already 
proposed by Ichimura et al. 2002 [4]. Besides the potato sequences 
from the PGSC Browser (Sotub and PGSC) the tree also includes 
four MAPKs from cv. Rywal: StMAPK4_l, StMAPK4_2, 
StMAPKO and StMAPK13. The sequences from Arabidopsis are 
AtMAPK3 (AT3G45640.1), AtMAPK4 (AT4G0 1370.1), At- 
MAPK5 (AT4G 11 330.1), AtMAPK6 (AT2G43790. 1), At- 
MAPK10 (AT3G59790.1), AtMAPKl 1 (AT1G01560.1), At- 
MAPK12 (AT2G46070.1) and AtMAPKl 3 (AT1G07880.2). 
The scale bar indicates the branch length that corresponds to 
0.06 substitutions per site. 
(TIF) 

Figure S3 Expression pattern of StMKK6 and its 
confirmed targets StMAPK4_2, StMAPK6 and StMAPK13 
in mock-inoculated plants. Log2 of normalized signals for 
mock-inoculated plants 1 day post inoculation are shown. The 
expression is shown for four mock treated Rywal (R) and NahG- 
Rywal (nah) plants. In both sets of plants the expression of the 
MKK6 interacting MAPKs is higher than of SxMKK6. Calculated 
are differences in the expression between Rywal and NahG-Rywal 
plants. 
(TIF) 

File SI Table SI in File SI. Oligonucleotides used for 
fragment isolation, qPCR and cloning. In primers for 
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cloning, the underlined parts of sequences are complementary to 
the destination plasmids. Table S2 in File SI. Validation of 
StMKK6 microarray (jiarray) results by real-time PCR 
(qPCR). To validate the microarray results for StMKK6, its 
expression was analysed by qPCR in the same RNA samples as for 
the microarray analysis. Log2 of ratio between virus- and mock- 
inoculated plants in cv. Rywal and NahG-Rywal 1, 3 and 6 days 
after PVY inoculation for StMKK6 obtained by both methods are 
shown. Statistically significant values (p<0.05) are marked with 
bold. Table S3 in File SI. Predicted regulatory domains of 
StMKKS native promoter from S. tuberosum, cv. Rywal. 
The promoter sequence is 899 bp-long and analysed with 
PlantCare software [36]. Table S4 in File SI. Predicted 
regulatory domains of StMKK6 native promoter from S. 
tuberosum, cv. Sante. The promoter sequence is 247 bp-long 
and analysed with PlantCare software [36]. Table S5 in File SI. 
Expression values of StMKK6 in control and SA-treated 
potatoes cv. Rywal and NahG-Rywal. Two biological 
replicates per treatment were analysed. Relative expression values 
and fold-changes (compared to Rywal control) are shown in the 
table. Differences between control and SA-treated plants were 
statistically evaluated by t-test. Table S6 in File SI. Compar- 
ison of StMKKS basal expression values in Rywal and 
NahG-Rywal plants. Two biological replicates per treatment 
were analysed. Relative expression values and fold-changes 
(compared to Rywal control) are shown in the table. Differences 
between Rywal and NahG-Rywal plants were statistically 
evaluated by t-test. Table S7 in File SI. Subcellular 
localisation prediction of StMKK6 and AtMKK6. Amino 
acid sequence of AtMKK6 (GenBank accession number 
NM_125041.2) and StMKK6 (GenBank accession number 
KF837 129.1) proteins were used as query in PredictProtein 
service. The localisation for each was predicted by three different 
prediction algorithms: PROSITE, LOCkey and LOCtree. For 
each subcellular prediction a confidence to the prediction is given. 
(DOC) 

File S2 Gene expression of potato MKK family in 
different tissues and after several stress treatments. 

For each MKK, the RNA expression signals (FPKM values), based 
on the RNA sequencing of double monoploid Solanum tuberosum 
Group Phureja clone DM 1-3 (DM) [40] and Genome sequence 
and analysis of the tuber crop potato [30] were collected in the 
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